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?H Abstract 

0\ We introduce results of computer experiments on information processing in 

a hexagonal array of vesicles filled with Belousov-Zhabotinsky (BZ) solution 
in a sub-excitable mode. We represent values of Boolean variables by ex- 
citation wave-fragments and implement basic logical gates by colliding the 
wave-fragments. We show that a vesicle filled with BZ mixture can imple- 
ment a range of basic logical functions. We cascade BZ-vesicle logical gates 
into arithmetic circuits implementing addition of two one-bit binary num- 
^ bers. We envisage that our theoretical results will be applied in chemical 

^ laboratory designs of massive-parallel computers based on fine-grained com- 

O partmentalisation of excitable chemical systems. 

0^ 
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1. Introduction 

X 

A reaction-diffusion computer is a spatially extended chemical system 
^ which processes information using interacting growing patterns of excitable 

and diffusive waves |3|. In react ion- diffusion processors, both the data and 
the results of the computation are encoded as concentration profiles of the 
reagents. The computation is performed via the spreading and interaction 
of wave fronts. All reaction-diffusion computers can be roughly split into 
two types: geometrically-constrained and free-space [6]. In geometrically- 
constrained computers excitation or chemical waves propagate in the 'hard- 
ware' of substrate defined channels and interact at the junctions between 



Preprint submitted to Elsevier 



June 11, 2010 



the channels. A great deal of experimental reaction-diffusion computing cir- 
cuits were implemented in a geometrically-constrained media: logical gates 
for Boolean and multiple- valued logic [261 ESI HZl ES], many- input logical 
gates [m [16], counters [13], coincidence detector [12], detectors of direc- 
tion and distance to a source of periodic oscillations [lU EO] and inductive 
memory [23] . 

Free-space reaction-diffusion computer^ do not have any underlying ar- 
chitecture. A computing medium is uniform and homogeneous, the only 
externally evoked disturbances of the medium's characteristics implement 
computation. In majority of cases computation in free-space computers is ex- 
ecuted using principles of a collision-based computation [^ . A collision-based 
computer employs mobile compact finite patterns and mobile self-localized 
excitations to represent quanta of information in active non-linear mediums. 
Information values, e.g. truth values of logical variables, are given by either 
absence or presence of the localizations or other parameters of the localiza- 
tions. The localizations travel in space and when collisions occur the result 
can be interpreted as computation. There are no predetermined stationary 
signal channels (wires), a trajectory of the travelling pattern is a momen- 
tarily wire. Almost any part of the medium space can be used as a wire. 
Localizations can collide anywhere within a space sample, there are no fixed 
positions at which specific operations occur, nor location specified gates with 
fixed operations. The localizations undergo transformations, form bound 
states, annihilate or fuse when they interact with other mobile patterns. In- 
formation values of localizations are transformed as a result of collision [2]. 

So far no advanced arithmetical circuits have been implemented in free- 
space reaction-diffusion chemical processors. Some preparatory steps are 
done however. They include simulation and experimental laboratory real- 
isation of basic logical gates [H HI [27] and generators of mobile localiza- 
tions, excitation wave-fragments (they can play a role of constant True) 
in light-sensitive Belousov-Zhabotinsky (BZ) medium [8] and adaptive de- 
sign of simple logical gates using machine learning techniques [2H]- The 
main reason for such slow progress is instability of wave-fragments that 
either collapse or expand without external control. One-bit half-adder in 
geometrically-constrained light-sensitive Belousov-Zhabotinsky medium was 
simulated in [5], however the implementation required dynamical update of 
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illumination level to prevent wave-fragments from collapsing or expanding. 
In this present paper we try to overcome the problem of excitation wave- 
fragments instability via combining geometrically-constraining and collision- 
based approaches. 

The paper is structured as follows. Two- variable Oregonator model of BZ 
medium is introduced in Sect. [2j Section |3] outlines principles of compart- 
mentalisation of BZ medium in planar discs and packing of the discs into a 
two-dimensional hexagonal lattice. Typology of interactions between wave- 
fragments is presented in Sect.|4j In Sect.[5]we introduce a cellular-automaton 
model of interacting wave-fragments. We show what types of logical gates 
can be realised via collision of wave-fragments in Sect. [6] These gates are 
employed to construct a one-bit full adder in Sect. [7| Future developments 
of the approach are outlined in Sect. [8j 



2. Oregonator model 

We use two- variable Oregonator equation [TO] adapted to a light-sensitive 
Belousov-Zhabotinsky (BZ) reaction with applied illumination |7]. 




In framework of BZ reaction the variables u and v represent local concen- 
trations of activator, or excitatory component, and inhibitor, or refractory 
component. Parameter e sets up a ratio of time scale of variables u and f , 
g is a scaling parameter depending on rates of activation/propagation and 
inhibition, / is a stoichiometric coefficient. Constant is a rate of inhibitor 
production. In light-sensitive BZ represents rate of inhibitor production 
proportional to intensity of illumination ([T]). 

To integrate the system we use Euler method with five-node Laplace 
operator, time step At = 0.001 and grid point spacing Ax = 0.25, e = 0.0243, 
/ = 1.4, q = 0.002. In some cases it was enough to speed up simulation 
by increasing At to 0.005 without loss of phenomenological accuracy. The 
parameter characterizes excitability of the simulated medium: the medium 
is excitable and exhibits 'classical' target waves when = 0.05 and the 
medium is sub-excitable with propagating localizations, or wave-fragments. 
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when = 0.0766. Time lapse snapshots provided in the paper were recorded 
at every 150 time steps, we display sites with u > 0.04. 

3. Compartmentalisation 




(c) 

Figure 1: Time lapsed snapshots of wave-fragments propagating in simulated BZ medium 
in (a) near lower threshold of excitability, (p — 0.079, (b) sub-excitable mode, (j> ~ 0.07905, 
(c) non-excitable mode, (f> = 0.0791. The media were perturbed by rectangular north-south 
elongated domains of excitation, 3 x 40 sites in state u = 1.0. Sites of initial segment-wise 
perturbation are shown by arrows. Grid size is 1125x250 nodes. Time lapsed snapshots 
are recorded every 150*'* step of numerical integration. 

Given initial asymmetric excitation domain a wave-fragment is formed. 
The fragment's velocity vector is a normal to longest side of the pertur- 
bation domain. In Fig. [T] a medium is excited by rectangular domains of 
perturbed sites. The perturbation domains are elongated along north-south 
axis therefore wave- fragments generated propagate west and east. Depending 
on medium's excitability wave- fragments may expand (Fig. [l^), keep their 
shape for a long time (Fig. ^) or collapse (Fig. 
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In ideal situation, assuming that wave-fragments keep their shapes indef- 
initely, we can implement a collision-based computing circuit of any depth 
subject to space availability. However in reality, particularly in conditions 
of chemical laboratory experiments, wave-fragments are very unstable. It is 
almost impossible to keep a medium at the precise level of sub-excitability 
and almost any wave-fragment will expand (Fig. [l^) or collapse 
The only solution, for controlling shape of wave-fragments, sug- 
gested so far was to change excitability of a medium periodically, as achieved 
in [25j in experiments with light-sensitive BZ medium. When a wave-fragment 
expands experimenters increase illumination level thus decreasing the medium's 
excitability. When the wave-fragment starts to collapse the illumination is 
decreased, the medium's excitability increases and the wave-fragment ex- 
pands. This approach does not look feasible for collision-based computing |2] 
with excitation wave-fragments, because it requires precise synchronisation 
of cycles of decreasing and increasing excitability. 

We can overcome the problem of wave-fragment instability by the subdi- 
viding computing substrate into interconnected compartments — BZ- vesicles 
— and allowing waves to collide one with another only inside the compart- 
ments. Each BZ-vesicle has a membrane impassable for excitation. A pore, 
or a channel, between two vesicles is formed when two vesicles come into 
direct contact. The pore is small such that when a wave passes through 
the pore there is not enough time for a wave to expand or collapse before 
interacting with waves entering through other pores, or sites of contact. 

A spherical compartment — BZ-vesicle — is the best natural choice. It 
also conforms to experimental results on encapsulating excitable medium in 
a lipid membrane and allows for effortless arrangement of the vesicles 

into a regular lattice. 

In this present paper we deal with a two-dimensional medium where vesi- 
cles are tightly packed into a hexagonal lattice (Fig. [2^). Each vesicle has 
radius 27 grid units, a pore has radius 5 units. The size of the vesicle and 
pore are critical in modulating the beam of the excitation waves. Changes 
in vesicle size will require further adjustment of parameter (Eq. [T]) in order 
to maintain the ideal wave spread within the vesicle. 

Each internal vesicle x has six closest neighbours xq - ■ ■ x^. Medium inside 
BZ-vesicles is sub-excitable while the medium outside the vesicles are non- 
excitable. Excitation wave-fragments propagate inside vesicles only. The 
waves pass from one vesicle through pores, sites of contacts. 
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(a) (b) 

Figure 2: Basic architecture element of BZ compartmentalized medium: (a) structure of 
disc-compartment neighbourhood, (b) excitation propagation on a regular arrangement of 
vesicles. 
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Proposition 1. BZ-vesicles packed in a hexagonal lattice provide a medium 
for directional transmission of information in a spatially extended non-linear 
excitable medium, where each vesicle acts as a signal amplifier and a pore, 
connecting two vesicles, acts as a focusing lens. 

We excite the medium inside one vesicle, say x, with a symmetric disc- 
shaped domain of perturbation centered in one of the vesicles. A circular 
excitation wave is formed. The circular wave reaches boundary of vesicle 
X, and enters pores between x and its neighbours Xq - ■ ■ x^ (Fig. |2|d). When 
excitation leaves a pore it has a form of a propagating wave-fragment. Due 
to the near-excitability of the medium (Fig. [l^) wave-fragments entering 
vesicles xq - ■ ■ x^ start to expand. However they do not have enough space to 
form a proper circular wave. Therefore they leave second-order neighbours 
of X via one exit pore per vesicle (Fig. |2]d). 

4. Interaction between wave-fragments 

Assumption 1. All constructs presented in the paper assume system of BZ- 
vesicles is fully synchronized. Waves enter any single vesicle simultaneously. 

This is a very strong and somewhat unrealistic assumption. However we 
use it as a starting point to build first reliable models of BZ- vesicle computers. 



Proposition 2. Let several wave-fragments enter a vesicle. If at least two 
wave- fragments have opposite velocity vectors all wave-fragments annihilate. 
Otherwise the wave-fragments merge and velocity vector of newly formed 
wave-fragment is a sum of velocity vectors of incoming wave-fragments . 

Constructive proof is illustrated by representative scenarios of wave in- 
teractions inside a single vesicle (vesicle x in Fig. [2^) are shown in Fig. |3} If 
just one neighbour of vesicle x is excited the vesicle x acts as a conductor 
and signal amplifier: the wave simply passes through the vesicle x slightly 
increasing in size and exits through the pore opposite to the wave's entry 
pore. Thus, in Fig. |3^ north-west neighbours (vesicle xq in Fig. [2^) is ac- 
tivated. Excitation enters vesicle x, wave-fragment travelling south-east is 
formed (Fig. |3^). The wave-fragment is transmitted to south-east neighbour 
(vesicle x^ in Fig. [2^) of vesicle x. 

There are three possible scenarios when two neighbours of vesicle x are ex- 
cited. In a situation when excited neighbours of x are also each others closest 



7 




8 



neighbours (e.g. vesicles Xq and Xi in Fig. [2^) the wave-fragments generated 
by them merge inside vesicle x (Fig. [3|3). The velocity vector of a newly 
formed wave-fragment is a sum of vectors of two original wave-fragments. 
Vectors of original wave-fronts orientate towards exit pores opposite to the 
excitation entry pores. The vector of a newly formed wave-fragment aims 
between the pores, hence no excitation leaves the vesicle x. For example, 
in (Fig. [3]d) north-west and north-east neighbours of vesicle x are excited. 
Two wave-fragments enter vesicle x: one fragment travels south-east, another 
south-west. The wave-fragments merge and form a new wave-fragment which 
travels south. This fragment collides with the part of vesicle x's wall lying be- 
tween south-west and south-east pores. The fragment becomes extinguished 
in result of the collision. 

If excited neighbours of vesicle x are separated by another neighbour of 
X (e.g. vesicles xq and X2 in Fig. [2^) then vector of the newly formed wave- 
fragment (result of merging two input wave-fragments) points exactly to one 
exit pore (e.g. pore connecting vesicle x and vesicle 0:4 in Fig. [2^). In example 
(Fig. [3]:) north-west and east neighbours of vesicle x are excited. Two wave- 
fragments enter vesicle x: one travels south-east another travels west. They 
merge and form a wave-fragment travelling south-west. This fragment leaves 
vesicle x for its southwest neighbour (vesicle X4). 

If wave-fragments travel towards each other (Fig. |3]i) they collide and 
annihilate in result of the collision. Thus no excitation leaves vesicle x. 
Scenarios of collision between three wave-fragments are shown in (Fig. [3^-g). 

In scenarios illustrated in (Fig. [3}i-k) at least two wave-fragments un- 
dergo head-on collision and all input wave-fragment annihilate. The situa- 
tion (Fig. [sji) is similar to (Fig. [sj^) with the only difference that three not 
two waves collide, however the resulting wave-fragment is the same: it travels 
south-west and excites south-west neighbour (vesicle X4) of vesicle x. 

Configuration of excited neighbours shown in Fig. |3]F results in annihila- 
tion of all three incoming wave-fragments. Three wave-fragments entering 
vesicle x travel south-east, south-west and north-west. Wave- fragment trav- 
elling south-east collides head-on with wave- fragment travelling north-west. 
Both wave-fragments annihilate as a result of the collision. At the same 
wave-fragment travelling south-west collides in both south-east and south- 
west fragments. The south-west travelling fragment also annihilates (Fig.jsjf). 

Situation shown in (Fig. [3^) does correspond to direct head-on collision 
between wave- fragments. Three incoming wave-fragments travel south-east, 
west and north-east. A sum of velocity vectors of these wave- fragments is 
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nil, therefore all fragments annihilate and no wave-fragments leave vesicle x. 

In the remaining scenarios (Fig.[3}i-k) at least two incoming wave-fragments 
experience head-on collisions and all wave-fragments annihilate. 

5. Cellular automaton representation 

We simulate interactions of wave-fragments in BZ-vesicles using two- 
dimensional cellular automaton. The cellular automaton is a hexagonal array 
L of cells. A cell is a finite state machine which updates its states in discrete 
time depending on states of its six closest neighbours. All cells update their 
states simultaneously. 

A cell takes eight states from set Q = {o, — , \, \, — )■}. They 

are :- resting state (o), when the cell is ready to be excited, refractory state 
(— ), when the cell does not react to states of its neighbours, and six excited 
states representing wave-fragments. The states representing wave-fragments 
are coded as \ (south-east travelling wave), (south-west travelling wave), 
^ (westward travelling wave), \ (north-west travelling wave), (north- 
east travelling wave) and — >■ (eastward travelling wave). The excited state of 
a cell indicates what type of excitation wave-fragment is leaving the cell. 

Let s : L X Q — )• {0, 1} defined as follows. For every x with neighbourhood 
u{x) = {xo, xi, . . . , X5} (Fig. [2^): s(xo)* = 1 if Xq =\ and s(xo)* = 
otherwise; si^XiY = 1 if =^ and s(xi)* = otherwise; ...; 5(^5)* = 1 
if x\ =— )■ and s{x^y = otherwise. Cell-state transition rule is defined as 
follows: 
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Example of the cellular automaton evolution from initial random config- 
uration is shown in Fig. |4j In case of absorbing boundary conditions any 
random initial configurations leads to 'empty' global configuration where all 
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(c) t = 8 (d) i = 9 

Figure 4: Configurations of tlie cellular automaton from initially random configuration: a 
cell gets on of six wave- fragment-states with probability 0.1 and gets a resting state oth- 
erwise. Boundary are absorbing. Resting cell-states are shown by small circles, refractory 
states are not shown. Cell-states representing wave-fragments are shown by large circles 
with vectors pointed south-east, south-west, west, north-west, north-east, east. 
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Figure 5: Traces of seven steps of cellular automaton development with eleven initial sub- 
configurations corresponding to types of wave collisions in Fig. [3j Cell-states representing 
wave-fragments are shown by large circles with vectors pointed south-east, south-west, 
west, north-west, north-east, east. To show the trace of the cellular automaton develop- 
ment we do not execute transition to refractory state but leave a cell in its wave-fragment 
state instead. Such approach is appropriate when signals trajectories do not cross each 
other. 

cells are in resting states. This is because wave-fragment either annihilate 
each other (in majority of local configurations) or produce less new wave- 
fragments than those involved in a collision. 
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Figure 6: Two and three input Boolean gates implemented in BZ-vesicles and the cellular 
automaton model. Solid lines indicate inputs, dashed lines - outputs. 



Examples of collision scenarios of Fig. [3] represented in cellular automaton 
configurations are shown in Fig. [5j The cellular automaton can be used as a 
'rapid prototyping' tool to design logical circuits in BZ-vesicle arrays. 

6. Boolean gates 

We adopt principles of collision-based computing [2] to design logical 
gates in BZ-vesicle arrays. We assume travelling wave-fragments represent 
values of Boolean variables: presence of a wave-fragment in a specified site 
corresponds to logical True, absence to logical False. When two wave- 
fragments collide they annihilate or merge into a new wave-fragment. Newly 
formed wave-fragments represent a logical conjunction. A wave-fragment 
travelling along its original (undisturbed) trajectory gives us a conjunction 
of a logical variable represented by this wave-fragment with negation of logical 
variables represented by another wave-fragment. 

With regards to usages of space by incoming and outgoing wave-fragments 
collision-based gates can be classified on gates where input and output tra- 
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jectories do not overlap (Fig. |6pLbd) and gates where we may record output 
signals from inputs (Fig. |6^bdf). 



6.1. Gates with separate input and outputs 

Gate (x, y) — )■ {xy, xy) (Fig. |6^) is implemented in the situation when 
wave-fragments incoming to a vesicle originate from the vesicles which are 
closest neighbours to each other (Fig. [3|d and Fig. [sjo). South-east travelling 
wave-fragment represents x and south-west travelling wave-fragment repre- 
sents y. If just one wave-fragment, say is present then it continues its 
course undisturbed, thus its outgoing trajectory symbolises xy. When two 
wave-fragments are present (x =true and y =true) they collide, form a 
new wave-fragment which annihilate by colliding in a vesicle wall between 
pores. 

In collision scenario (Fig. [sj: and Fig. [s]:) wave-fragments x and y do 
not annihilate but form a new wave-fragment which travels along its own 
trajectory. This new wave-fragment represents conjunction xy. Thus we 
have a two- input three-output logical gate (x,?/) — {xy, xy,xy) (Fig. |6|d). 

Arrangement of colliding wave-fragments shown in Fig. |3^ and Fig. |5^ 
implements three-input three-output gate {x, y, z) — )• {xyz, xyz + xz, xyz) 
(Fig. [6]i). By assigning constant True to one of the inputs we can trans- 
form the gate to negation of exclusive disjunction. An implementation of the 
gate (x, y) — )■ (x © y) with collision of wave-fragments is shown in Fig. [7| 
This example also demonstrates that by exciting input vesicles with omnidi- 
rectional waves we multiply original signals, and copies of the signals can be 
employed in further gates. 

6.2. Gates with spatially overlapping inputs and outputs 

Gate shown in Fig. ^ is the simplest amongst gates with spatially over- 
lapping trajectories of input and output wave- fragments. If ?/=False then no 
wave-fragment is present in input y. Therefore wave-fragment representing 
X =True travels along its original trajectory undisturbed and eventually 
appears on the output xy. The three-input, three-output gate in Fig. |6^ 
produces its result xyz, xyz, xyz on the output trajectories co-aligned with 
input trajectories for signals x, y and z. 

Gate Fig. [6]F has three inputs and six outputs. Three outputs — xyz, xyz, 
xyz — are represented by wave-fragments travelling outward the collision site 
but along the input trajectories of signals x, z and y. Three other outputs 
have their own trajectories — xyz, xyz, and xyz. 
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y y 




Figure 7: Example of gate x ®y. The figure shows time lapse images of propagating 
excitation wave-fragments: (a) experimental setup with omnidirectional excitation of three 
sites of BZ-vesicle array, (b) selected local configuration of the gate. Constant True is 
shown by '1'. 
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While it may be convenient to employ gates with spatially overlapping 
inputs and outputs, e.g. when space is an issue, a proper functioning of the 
gates will require more advanced temporal coordination of signals. 



7. Towards binary adder 

We employ some of the gates with spatially separated input and output 
trajectories to design components of one-bit binary adder. The example 
provided is for illustrative purposes: we did not optimize structure of logical 
circuits presented, neither we were concerned with number of instances of 
input variables used. 

A binary adder based on BZ-vesicles uses collisions between propagating 
wave-fragments to perform addition of three one-bit binary numbers x, y 
and z. The signal z depict Cm- There are many versions of particular 
implementations of full one-bit adder. We adopted the most common one. 



where the adder outputs sum of signals S = {x (B y) (B z (Fig. [8]& Fig. 10 1) 
and carry out value Cout = xy + z{x © y) (Fig. [9]& Fig[Io|3). 

Positions of inputs in Figs. |8]and|9]are shown by thick solid circles. If an 
input equals False the circle contains only central dot, if the input equals 
True the circle contains a vector indicate the sites state. Result of sum 
circuit is represented by wave-fragment travelling south-east while output 
of carry out circuit is represented by wave- fragment travelling east. The 
circuit calculating S is packed in array of 6 x 10 BZ-vesicles, and the circuit 
calculating Cout occupies as sub-array of 7 x 13 BZ-vesicles. Let us discuss 
functioning of circuit S (Fig. Isl) in details. 



7. 1 . Sum 

If all three inputs are False no waves are initiated in sites and 
and W7 (Fig. [8^ and 10 1). The wave- fragments representing constant True 



are placed in sites wi and W2 (south-east travelling wave-fragments), W4 and 
wq (wave- fragments travelling east). Wave- fragment W2 collides with wave- 
fragment W4, both wave- fragments annihilate (step three of cellular automa- 
ton simulation. Fig. |8^). Wave- fragments Wi and Wq travel a bit further but 
collide one with another and annihilate at fifth step of automaton simulation. 
All outputs are False therefore. 

If carry in value z is True and other input values are False (Fig. |8^ and 



10 d) then north-east traveling wave annihilates east travelling wave w^. 
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(d) X = 0,2/ = 1,2 = 1 (e) X = 1,2/ = 0,z = (f) x = 1, 2/ = 0, z = 1 
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(g) a; = l,y = 1,^ = (h) X = 1,2/ = l,z = 1 

Figure 8: Calculation of sum S — {x(Sy)(Bz,z = Ci„, in a cellular automaton model of 
BZ-vesicle hexagonal array: (a)-(h) traces of cell-states, representing wave-fragments, for 
all possible combinations of input values (Fig. 10). Input cells are shown by circles with 
thick lines. 
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(g) x=^l,y = l,z = (h) X = 1, y = 1, z = 1 

Figure 9: Calculation of carry out value Cout — xy + z{x ® y), z ~ Ci„ in cellular au- 
tomaton model of BZ-vesicle hexagonal array: (a)-(h) traces of cell-states, representing 
wave- fragments, for all possible combinations of input values (Fig. 10). Input cells are 
shown by circles with thick lines. Eastward travelling wave- fragments in 3rd and 13th 
rows in (ef) are unused byproducts of the circuit. 
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Figure 10: Scheme of the circuit, (a) Sum (Fig.ls]) & (b) Carry out (Fig. [9]). 
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Therefore wave Wi continues travelling south-east, thus representing True 
value of (x © y) © z. 

In situation {x, y, z) = (False, True, False) wave collides with 
wave W4 (Fig. [sj: and [lO^) , both wave-fragments annihilate. Therefore wave- 
fragment W2 continues undisturbed its travel to south-east where it collides 
with wave-fragment Wq. With both waves representing east travelling con- 
stant True cancelled wave- fragment wi reaches an output. Similarly for 
input {x,y,z) = (True, False, False) wave-fragments and W4, and W2 
and wq annihilate in collision with each other (Fig. [8^ and [lO^) . Therefore 
wave-fragment Wi travels undisturbed. Other combinations of inputs can be 
considered in the similar manner. 

7.2. Carry out 

Computation starts in a group of sites marked Wi, . . .w^m (Fig.[lO|3): site 
W2 represents x, site represents y, and sites Wi and represent constant 
True. The sub-circuits' output is x © y. Wave- fragments Wi, travelling 
south-east, and W3, travelling east, are always present in the system. If only 
one of inputs x or y has True value, e.g. south-east travelling wave-front W2 
in Fig. [9^f), then wave-fragment representing this input collides with wave- 
fragment W3 (constant True) and both wave-fragments annihilate. If both 
inputs X and y are True (Fig. [9^h) then waves W2 and W4 collide with each 
other and wave 1^3, merge altogether and produce new wave travelling east. 
This new wave (seen in node (6, 5) in (Fig. [9^h)) collides with wave wi, both 
wave- fragments annihilate. Thus, sub-circuit {wi, . . . , W4) computes x ®y. 

Output wave-fragment of sub-circuit x ® y collides with wave-fragment 
W5, which represents carry in value z. These wave-fragments collide at angle 
120°, therefore they merge and produce new wave-fragment (travelling east) 
when collide. See examples in Fig. |9]if. Thus intermediary result z{x®y) is 
calculated. 

Two small sub-circuits — Wq and w-j, and and - are arranged 
symmetrically north and 6,6south of trajectory of wave-fragment, which 
represents z{x (By). Each of the sub-circuits (^6,^7) and (^§,^9) pro- 
duces xy. The final wave-fragment wiO is only produced when either only 
wave- fragment z{x (B y) travels east (Fig. [ojif), or all three wave-fragments 
— wave-fragment xy (output of sub-circuits (^6,^7)) travelling south-east, 
wave-fragment xy (output of sub-circuits {ws,Wg)) travelling north-east, and 
wave- fragment z{x © y) collide (Fig. |9|i). 
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8. Summary 

Belousov-Zhabotinsky (BZ) medium in a sub-excitable mode exhibit the 
rich spatio-temporal dynamics of mobile self-localisations, wave-fragments 
behaving as distant analogs of dissipative solitons. In an ideal situation 
a wave-fragment keeps its shape and travels in a predetermined direction 
until a collision with another wave-fragment. Therefore wave- fragments can 
be used as a quanta of information, or analogs of billiard vesicles [11], in 
developing collision-based computing circuits [2] . Our attempts to implement 
large scale logical circuits in experimental chemical laboratory conditions 
show that wave-fragments are very unstable and do not conserve their shape 
for a long time without external control [271 |H]. 

To overcome the problem of wave-fragment instability we decided to 
achieve a compromise between geometrically-constrained and free space ap- 
proaches in designing reaction-diffusion chemical computers. In computer 
simulation experiments we designed a lattice of two-dimensional vesicles (BZ- 
vesicles) filled with BZ reaction solution. The BZ-vesicles are in close contact 
with each other, via tiny pores, and thus excitation wave can pass form one 
vesicle to its closest neighbour. Wave-fragments keep their shape or just 
slightly expand inside each BZ-vesicle. When a wave-fragment passes from 
one vesicle to another it contracts, due to the restricted size of the connecting 
pore. 

In computer experiments we demonstrated that each BZ-vesicle, in a 
hexagonal arrangement of vesicles, can act as a logical gate with at most six 
inputs and six outputs. All basic operations of a Boolean logic are imple- 
mented via collision between wave-fragments in a vesicle and thus a lattice 
of BZ-vesicles is a computationally universal substrate. To show that BZ- 
vesicles pose not only a purely theoretical interest we cascaded the logical 
gates, mostly gates with two inputs and one or two outputs, in logical circuits 
which implement addition of two one-bit binary numbers. Functioning of the 
circuits was verified in cellular-automaton models of BZ-vesicle array. 

Implementation of theoretical designs in chemical laboratory experiments 
would be the next step of our studies in fine-grained compartmentalised ex- 
citable chemical systems. Some experimental works have been done, Kita- 
hata and colleagues provided amazing experimental evidences of excitation 
wave-travelling inside a tiny sphere and even movement of BZ-vesicle partly 
controlled by wave propagation [211 122] • Also, BZ reaction mixture can be 
encapsulated via dispersion of the mixture in a water-in-oil micro-emulsion 
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with surfactant: BZ reagents becomes enclosed in a mono-layer of anionic 
surfactant. 'Unlimited' energy supply from solution surrounding BZ-vesicles 
enables reactivity of BZ droplets to be sustained for extended periods of 
time [9]. 
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